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INTERACTION  BETWEEN  A  SUPERSONIC  STREAM  AND  A  PARALLEL 
SUBSONIC  STREAM  BOUNDED  BY  FLUID  AT  REST 
By  Herbert  S.  Ribner  and  E.  Leonard  Arnoff  • 


SUMMARY 

In  a  simplified  inviscid  model  of  shock-wave  boundary-layer  inter¬ 
action,  Tsien  and  Finston  have'  replaced  the  boundary  layer  by  a  uniform 
subsonic  stream  bounded  on  one  side  by  a  solid  wall  and  on  the  other 
side  by  the  interface  with  a  uniform  supersonic  stream  of  semi-infinite 
extent.  Among  other  things,  this  model  fails  to  simulate  the  separated 
region  or  "dead-air"  bubble  that  generally  appears  in  a  laminar  boundary 
layer  subjected  to  an  oblique  incident  shock  wave  of  moderate  strength, 
pi  order  to  introduce  a  main  feature  of  such  a  dead-air  region,  the  model 
has  been  modified  herein  by  replacing  the  solid  wall  by  an  interface  with 
fluid  at  rest. 

The  presence  of  the  boundary  layer  sandwiched  between  the  outer 
supersonic  flow  and  the  dead-air  region is  found  scarcely  to  modify  the 
shape,  in  the  vicinity  of  the  shock,  of  the  expansive  "corner"  turn  that 
would  exist  if  the  shock  were  incident  directly  on  the  dead-air  region 
without  the  intermediary  of  the  boundary  layer;  there  are  local  distor¬ 
tions  top  and  bottom,  but  these  are  reduced  to  negligible  amounts  several 
boundary-layer  thicknesses  to  the  left  or  right  of  the  effective  corner. 

In  support  of  a  phase  of  the  work  of  Lester  Lees,  it  is  concluded 
that  in  a  more  accurate  treatment  of  the  complete  region  of  shock 
boundary-layer  interaction  the  Prandtl  boundary -layer  equations  may  be 
applied  to  the  entire  extent  of  the  disturbed  boundary  layer,  applying 
as  a  boundary  condition  a  sudden  turn  of  the  displacement  surf, ace  through 
twice  the  shock  deflection  angle  at  the  point  of  shock  incidence.  There¬ 
by  the  flow  details  in  the  immediate  vicinity  of  the  shock  will  be  some¬ 
what  in  error,  but  the  over-all  features  of  the  interaction  are  capable 
of  being  given  correctly.  Present  unknown  elements  in  such  an  application 
appear  to  be  the  point  of  transition  from  laminar  to  turbulent  flow  and 
the  form  of  the  turbulent  equations  where  a  separated  bubble  exists.  < 
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INTRODUCTION 

The  complex  problem  of  the  interaction  between  shock  waves  and 
boundary  layers  in  supersonic  flow  is  not  yet  clearly  understood, 
despite  a  number  of  theoretical  studies  (references  1  to  7).  Among 
these,  the  analysis  of  Tsien  and  Finston  (reference  2)  is  singled 
out  here  for  further  examination.  In  the  Tsien-Finston  model  the 
boundary  layer  is  simulated  by  a  uniform  subsonic  stream  of  finite 
width  bounded  on  one  side  by  a  solid  wall  and  on  the  other  side  by 
the  interface  with' a  uniform  supersonic  stream  of  semi-infinite 
extent.  The  fluid  is  assumed  to  be  nonviscous  and  nonheat -conducting 
and  the  disturbances  are  assumed  small. 

The  boundary-layer  perturbation  calculated  in  reference  2  from  this 
model  has  the  general  character  shown  in  sketch  1.  The  experimental 
perturbation  for  a  laminar  boundary  layer  (e.g.,  references  8  and  9), 
on  the  other  hand,  has  the  general  character  shown  in  sketch  2 .  The 
resemblance  is  not  close.  A  striking  experimental  feature  is  the  pres¬ 
ence  of  a  separated  region,  which  appears  for  all  but  the  weakest  inci¬ 
dent  shock  strengths  (references  3  and  9).  The  separated  region  is  not  . 
simulated  in  the  Tsien-Finston  model  of  shock  boundary-layer  interaction, 
and  this  constitutes  an  important  shortcoming  of  the  model  for  laminar 
boundary  layers.  For  turbulent  boundary  layers  subjected  to  moderate 
shock  strengths,  on  the  other  hand,  there  will  be  no  separation,  and 
the  model  may  be  quite  good.- 

The  air  in  a  separated  bubble  is  substantially  at  rest.  An  oblique 
shock  incident  directly  on  a  dead-air  region,  without  the  intermediary 
of  a  boundary  layer,  is  known  to  reflect  as  an  expansion  wave  (see 
sketch  3) .  Note  that  the  interface  in  sketch  3  is  deflected  much  like 
the  potential  flow  just  above  the  actual  boundary  layer  (sketch  2),  a 
feature  completely  unduplicated  in  the  Tsien-Finston  model  (sketch  l). 

A  natural  inference  is  that  the  separated  region  dominates  the  behavior 
in  the  immediate  vicinity  of  the  shock,  and  the  upper  surface  behaves, 
in  that  vicinity,  substantially  as  if  the  entire  boundary  layer  were 
replaced  by  a  dead-air  region. 

This  tempting  inference  has  been  made  to  provide  a  boundary  condi¬ 
tion  in  the  work  of  Lees  (reference  6)  and  of  Ritter  (reference  l),  along 
with  other  simplifying  assumptions;  it  is  the  purpose  of  the  present 
paper  to  test  its  validity  on  a  simplified  model  applicable  only  to. the 
immediate  vicinity  of  the  shock.  For  this  purpose  a  modification  of 
the  Tsien-Finston  model  will  serve:  the  solid  lower  wall  of  sketch  1  is 
replaced  by  a  semi-infinite  region  of  air  at  rest,  or  "dead  air".  Thus 
the  boundary  condition  of  zero  vertical  perturbation  velocity  at  the 
wall  is  replaced  by  a  condition  of  constant  pressure  along  the  boundary 
of  the  dead-air  region;  in  all  other  respects  there  is  no  change  from 
the  original  model.  Experimentally  (reference  9),  an  approximation  to 
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the  condition  of  constant  pressure  is  found  over  a  substantial  portion 
of  the  separated  air  bubble  including  the  part  directly  under  the  inci¬ 
dent  shockj  the  ultimate  sharp  pressure  rise  generally  occurs  somewhat 
downstream  of  the  shock. 

It  has  been  implicit  in  the  discussion  that  attention  was  limited 
to  the  boundary  layer  adjacent  to  a  solid  body.  Also  of  interest,  how¬ 
ever,  is  the  boundary  layer  between  a  two-dimensional  supersonic  jet  and 
the  stationary  fluid  bounding  it.  The  present  modified  flow  model  will 
describe,  in  idealized  fashion,  the  interaction  of  an  oblique  shock  wave 
(such  as  may  emanate  from  a  nozzle  lip)  with  such  a  jet  boundary  layer. 

The  model  is,  in  fact,  better  adapted  to  such  a  fluid-to-fluid  boundary 
layer  than  it  is  to  the  fluid -to-wall  boundary  layer. 

The  mathematical  analysis  herein  of  the  modified  model  is  necessarily  > 
parallel  in  a  number  of  respects  to  the  original  work  of  Tsien  and  Finston. 
Nevertheless,  it  has  been  thought  worthwhile  to  provide  a  complete  account 
in  an  appendix  for  purposes  of  unity  and  clarity.  The  symbols  Used  in  this 
analysis  are  presented  in  appendix  A;  the  analysis,  in  appendix  B.  An 
abbreviated  version  of  the  analysis  is  given  in  the  main  text.  The  inves¬ 
tigation  was  conducted  at  the  NACA  Lewis  laboratory. 


SYNOPSIS  OF  ANALYSIS 

Flow  regions  and  governing  equations.  -  In  figure  1  are  shown  the 
supersonic  region  (region  lj  at  the  top,  the  subsonic  region  (region  2) 
in  the  middle,  and  the  dead-air  region  (region  3)  at  the  bottom.  A  weak 
shock  wave,  or  Mach  wave,  is  incident  on  the  interface  between  regions  1 
and  2.  The  governing  equations  for  the  perturbation  potential  are 

/  2\  d2<Pi  d2<P-> 

Region  1:  (l  -  1  — -g  .+  g  =  0,  Mx  >  1 

\  /  cbc  by^ 

/  2\  d2*2  d2?o 

Region  2:  ll  -  Mg  j  - g  +  g  ■  0,  Mg  <  1 

\  /  Sx  by 

In  region  1  the  perturbation  potential  includes  a  part  (incoming  wave) 
which  represents  directly  the  incident  Mach  wave  and  a  part  (outgoing 
waves)  which  arises  from  the  consequent  warping  of  the  interface  between 
regions  1  and  2. 

Boundary  conditions  at  each  interface.  -  The  boundary  conditions  are: 
(a)  The  perturbation  pressure  must  be  the  same  on  both  sides  of  an  inter¬ 
face  (taken  to  be  zero  at  the  lower  interface),  and  (b)  the  flow  inclina¬ 
tion  must  be  the  same  on  both  sides.  In  linearized- theory,  boundary  con¬ 
dition  (a)  is  essentially  a  condition  on  the  axial  perturbation  velocity 
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u,  and  boundary  condition  (b)  is  essentially  a  condition  on  the  trans¬ 
verse  perturbation  velocity  v .  At  the  lower  interface,  condition  (b) 
yields  no  information,  since  no  direction  can  be  assigned  to  the  zero 
velocity  in  the  bottom  region;  thus  (a)  and  (b),  as  applied  to  both 
interfaces,  yield  three  boundary  conditions. \ 

Solution  in  general  form.  -  The  general  solution  for  each  region  is 
expressed  in  the  form  of  a  Fourier  integral;  such  a  form  is  capable  of 
satisfying  arbitrary  boundary  conditions.  The  variable  of  integration 
is  A,  fl.nrl  terms  like  A-^(A)  cos  Ax,  B^(A)  sin  Ax,  and  so  forth,  appear; 

there  are  six  Fourier  coefficients  Ai( A),  B^(A),  E2(A),  F2(A),  ^(A), 
and  H2(A)  to  be  determined.  ' 

App1  "i cation  of  boundary  conditions.  -  The  velocity  components  in 
Fourier  integral  form  are  now  substituted  into  the  three  boundary  con¬ 
ditions;  equating  corresponding  coefficients  of  sine  terms  and  cosine 
terms  results,  in  effect,  in  six  equations.  Two  of  the  equations  yield 
E2(A)  =  F2(A)  =  0;  thus  there  remain  four  equations  which  are  solved 
simultaneously  for  the  four  remaining  Fourier  coefficients  A^A),  B1(A), 
G2(A),  and  H2(A). 

Solution  for  pressure  along  upper  interface.  -  The  perturbation 
pressure  along  the  upper  interface  is  proportional  to  u2  there;  the 
expression  is  of  the  form  of  a  Fourier  integral  involving  G2(A)  and 
H2(A),  which  are  now  known.  Evaluation  and  simplification  of  the 
integral  is  complicated  and  tedious;  the  final  result  is  expressed  in 
series  form  and  evaluated  numerically .  In  this  evaluation,  functions 
tabulated  by  Tsien  and  Finston  are  employed. 

Solution  for  displacement  of  each  interface.  -  The  local  displace¬ 
ment  t)(x)  of  an  interface  from  its  undisturbed  position  is  obtained  by 
integration  of  the  slope  dq/ dx;  this  slope  is  related  to  v2  by 

iasl2  ' 

dx  U2 

Like  u2,  v2  has  been  obtained  as  a  Fourier  integral,  involving  G2(X) 
and  H2(A),  and  G2  and  H2  are  now  known.  Again,  evaluation  and 
simplification  of  the  integrals  is  an  involved  process,  and  the  final 
displacements  are  exhibited  in  series  form  and  evaluated  numerically . 

RESULTS  AND  DISCUSSION 
Present  Modified  Tsien-Finston  Model 

Mach  number  parameter  9.  -  In  the  results  herein  the  Mach  numbers 
M^  send  Mg  of  the  respective  supersonic  and  subsonic  streams  appear  in 
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a  joint  parameter  0  as  well  as  individually..  For  a  given  value  of  0, 
M2  is  a  prescribed  function  of  Mj_  (see  equation  (B32)).  The  curves  of 
'  Mg  against  M]_  for  several  values  of  0  are  given  in  figure  2,  which 
has  been  reproduced  from  figure  3  of  reference  2.  Tsien  and  FinSton 
note  that  if  the  representative  Mach  numbers  for  the  boundary-layer, 
flow  and  the  free  stream  are  taken  to  be  0.8  and  2,  respectively,  then 
the  corresponding  value  of  0  in  the  flow  model  is  approximately  3jt/4. 

Pressure  distribution.  -  Curves  of  the  pressure  distribution  along 
the  upper  interface  calculated  from  equations  (B62)  and  (B63)  with  the 
aid  of  functions  tabulated  in  reference  2  are  given  in  figures  3(a), 

3(b),  and  3(c)  for  values  Jt/4,  jt/2,  and  3n/4,  respectively,  for  the 
Mach  number  parameter  0.  These  curves  exhibit  the  same  infinite  pres¬ 
sure  peak  at  the  point  of  shock  incidence  (£  =  0)  found  by  Tsien  and 
Finston  with  the  original  model.  A  similar  peak  was  found  by  Lighthill 
(reference  5)  in  his  model  with  continuous,  rather  than  stepwise,  veloc¬ 
ity  profile;  in  that  case  the  peak  occurred  inside  the  simulated  boundary 
layer  at  the  point  of  reflection  of  the  shock  from  the  sonic  line.  Such 
a  peak  will  be  modified  in  a  real  flow  and  probably  will  appear  as  a 
localized  region  of  high  pressure. 


Displacement  of  upper  and  lower  interfaces.  -  The  respective  dis¬ 
placements  of  the  upper  and  lower  interfaces  from  their  undisturbed 
positions  have  been  designated  Tjy  and  T)^;  both  are  functions  of  the 

longitudinal  distance  x.  Convenient  corresponding  nondimens ional 


parameters  are 


jt  % '  a 
2b(32  2t’  2bp2  2t’ 


and 


JtX 

2b(32‘ 


Curves  of  the  displacement 


against  longitudinal  distance  in  terms  of  these  nondimens ional  parameters 
are  plotted  in  figures  4(a),  4(b),  and  4(c)  for  0  equal  to  Jt/4,  Jt/2, 
and  3it/4,  respectively.  The  infinite  slope  at  x  =  0,  associated  with 
the  pressure  peak  there,  is  so  localized  on  the  drawings  as  to  be 
scarcely  distinguishable.  . 


The  scales  are  such  that  the  vertical  and  horizontal  displacements  .. 
of  a  given  interface  are  to  the  same  geometric  scale  for  the  case 
e  =  l/8  radian;  for  this  value  of  e  the  separation  of  the  interfaces 
is  to  correct  scale  when,  in  addition,  p2  =  Jt/20  for  0  =  Jt/4  (Mj_  =  2.5, 

Mg  =  0.988);  P2  =  Jt/l0  for  0  =  jt/2  (Mp  =  2.7,  Mg  =  0.95);  and  p2  =  Jt/6 
for  0  =  3jt/4  (M-]_  =  3.1,  Mg  =  0.85).  Numerical  values  of  the  ordinates 
and  abscissas  are  also  tabulated  in  tables  I  and  II. 


Comparison  with  supersonic  flow  above  still  air,  no  boundary  layer.  - 
The  curves  are  now  to  be  compared  with  the  conjecture  given  in  the  intro¬ 
duction:  namely,  that  in  the  simplified  model  of  figure  1  the  upper 
interface  behaves,  in  the  vicinity  of  the  shock,  substantially  as  if  the 
entire  boundary  layer  (represented  by  region  2)  were  replaced  by  a  dead- 
air  region  (an  extension  of  region  3);  that  is,  that  the  upper  interface 
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deflects  through  an  angle  2e  just  as  the  dead-air  interface  does  in  >  ' 
sketch  3.  The  comparison  shows  that  the  curves  in  figure  4  approximate 
the  conjectured  behavior;  both  interfaces  asymptotically  approach  the 
expansive  turn  through  an  angle  2e.  - 

For  each  interface,  the  final  asymptote  does  not  intersect  the 
initial  asymptote  (which  is  the  5-axis)  at  the  point  of  shock  incidence, 
5=0;  instead  the  point  of  intersection  lies  somewhat  downstream  thereof. 
The  point  of  intersection  of  the  two  asymptotes  may  be  thought  of  as  the 
effective  "comer"  for  the  turn  through  26;  its  displacement  amounts 
from  one  to  several  "reduced"  boundary -layer  thicknesses  (reduced  thick¬ 
ness  =  bp2)>  depending  on  the  value  of  the  Mach  number  parameter  0. 

The  local  deviation  from  the  asymptotic  behavior,  that  is,  from  the 
rectilinear  corner  turn  of  sketch  3,  is  limited  essentially  to  several 
reduced  boundary-layer  thicknesses  upstream  and  downstream  of  the  effec¬ 
tive  corner. 

Further  Generalization  of  Tsien-Finston  Model  and  Limitations 

The  flow  model  of  figure  1  is  supposed  to  simulate  only  one  aspect 
of  shock  boundary-layer  interaction,  namely,  the  local  effect  of  the 
separated  region  in  the  vicinity  of  the  point  of  shock  incidence.  An 
obvious  further  generalization,  later  found  to  be  of  questionable  util¬ 
ity,  is  sketched  in  figure  5.  Here  the  dead-air  region  is  limited  to  a 
finite  length  to  approximate  the  "separation  bubble,"  and  on  either  side 
thereof  the  solid-wall  condition  used  by  Tsienand  Finston  is  reapplied. 

When  points  S  and  R  are  both  more  than  several  reduced  boundary- 
layer  thicknesses  from  I,  the  point  of  shock  incidence,  their  effects 
in  the  neighborhood  of  X  would  be  expected  to  be  small.  This  consider¬ 
ation  serves  to  justify  the  use,  for  the  neighborhood  of  I,  of  the  model 
of  figure  1  in  which  S  and  R  are  effectively  removed  to  minus  and  plus 
infinity,  respectively. 

There  is  another  mentionable  difference  between  the  models  of  fig¬ 
ures  1  and  5:  In  figure  1  the  perturbation  pressure  in  the  dead-air 
region  is  taken  to  be  zero  (for  simplicity),  whereas  a  positive  pertur-. 
bation  pressure  is  intended  in  figure  5.  The  pressure  rise  in  figure  5 
is  effected  by  means  of  a  weak  shock  or  Mach  wave  incident  just  above  S, 
resulting  in  the  general  upward  slope  of  the  boundary  layer  between  S 
and  I.  This  more  realistic  positive  pressure  perturbation  in  the  dead- 
air  region  may  be  introduced  intp  the  model  of  figure  1  by  imagining  the 
weak  shock  associated  with  point  S  to  occur  at  minus  infinity  to  the 
left.  Then  the  entire  flow  field  will  have  an  upward  deflection  (say  e-jJ 
corresponding  to  the  flow  deflection,  through  this  shock.  As  a  conse¬ 
quence,  the  interface  deflections  of  figure  4  will  have  superposed  a 
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deflection  r,  =  «jX  (or  ^  h.  =  '  I?  >*  thls  is  eubntahtially 

as  though  each  figure  were  rotated  counterclockwise  through  the  angle 

el  • 

As  a  generalization  from  a  single  aspect  to  a  more  complete  picture 
of  shock  houndary-layer  interaction,  the  model  of  figure  5  has  serious 
shortcomings.  The  model  by  itself  is  insufficient  to  determine  the  posi¬ 
tions  of  separation  and  reattachment  S  and  Rj  for  this  purpose  the 
neglected  viscosity  will  have  to  he  brought  into  the  picture.  Moreover, 
between  I  and  R  the  actual  boundary-layer  flow  will  probably  become 
turbulent,  and  the  neglected  turbulent  momentum  exchange  across  both  of 
the  idealized  interfaces  will  markedly  influence  the  shape  of  the  effec¬ 
tive  interfaces  and  the  location  of  R.  The  transition  to  turbulent 
flow  will  probably  be  decisive  in  determining  ‘R  unless  the  incident 
shock  is  relatively  weak,  since  otherwise  the  large  abrupt  pressure  rise 
upon  reattachment  would  be  expected  to  reseparate  a  laminar  layer  but 
not  a  turbulent  one.  For  these  reasons  it  has  not  seemed  worthwhile  to 
make  a  serious  effort  to  obtain  an  analytical  solution  for  this  model. 


Applicability  of  the  Assumptions  of  Boundary- Layer  Theory 

On  the  present  model  of  shock  boundary-layer  interaction,  the  region 
in  which  the  upper  and  lower  interfaces  show  substantial  curvatures  and 
deviations  from  parallelism  has  been  shown  to  be  limited  to  one  or  two 
reduced  boundary-layer  thicknesses  on  either  side  of  the  point  of  shock 
incidence)  a  similar  result  is  found  for  the  original  Tsien-Finston 
model.  These  results  imply  that  it  is  only  in  this  region  that  the  com¬ 
ponent  pressure  gradient  dp/dy  is  significant.  Elsewhere  the  primary 
assumption  of  the  Prandtl  boundary-layer  theory,  namely,  that  dp/dy 
may  be  neglected,  appears  to  be  applicable. 

The  situation  seems  to  be  as  follows:  Both  models  of  shock  boundary 
layer  interaction  completely  neglect  viscosity,  but  include  the  inertial 
effects  in  both  x-  and  y-directions,  although  in  idealized  fashion.  The 
boundary-layer  theory,  on  the  other  hand,  includes  viscosity  and  the 
inertia  effects  in  the  x -direction  alone,  the  latter  with  idealization. 
The  results  of  these  models,  despite  their  neglect  of  viscosity  and 
other  idealizations,  suffice  to  delimit  the  small  region  directly  under 
the  incident  shock  in  which  the  inertia  effects  in  the  y-direction  must 
be  included.  It  appears  to  be  a  proper  inference  that  the  Prandtl 
boundary -layer  equations,  laminar  or  turbulent  as  the  case  may  be,  and 
with  provision  for  interaction  with  the  outer  flow,  may  be  applied 
upstream  and  downstream  of  this  region,  except  possibly  at  the  point 
of  reattachment  of  the  separated  bubble .  The  separated  bubble  itself  , ' 
except  that  portion  directly  under  the  shock,  is  included  in  the 
region  of  applicability. 
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Furthermore ,  the  extension  of  the  houndary-layer  equations  into  the 
forbidden  region,  with  the  assumption  of  a  sudden  turn  of  the  displace¬ 
ment  surface  through  an  angle  2c  at  the  shock,  appears  to  be  justified 
in  a  practical  analysis  where  accuracy  in  local  flow  details . at  the 
shock  is  secondary  to  the  over-all  picture;  this  conclusion  is  in  sup¬ 
port  of  a  phase  of  the  treatments  of  references  6  and  7 .  ^ If  the 

incident  shock  is  too  weak  to  cause  a  separated  bubble,  the  2e  turn 
at  the  shock  must  be  replaced  by  another'  condition. )  Present  unknown 
elements  in  such  an  application  appear  to  be  the  choice  of  the  point  of 
transition  to  turbulent  flow,  if  the  boundary  layer  is  not . initially tur¬ 
bulent,  and  equations  relating  pressure  distribution  and  displacement 
thickness  for  a  separated  region  of  a  turbulent  or  transitional  boundary 

layer. 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  September  19,  1952. 


TNote  added  in  proof:  In  a  newly  published  paper  ("A  Mixing  Theopr  for  _ 
the  Interaction  Between  Dissipative  Flows  and  Nearly  Isentropic  Streams, 
Jour.  Aero.  Sci.,  vol.  19,  no.  10,  Oct.  1952,  pp.  649-676)  Crocco  and 
Lees  have,  in  effect,  supplied  the  missing  relation -between  pressure  dis¬ 
tribution  and  displacement  thickness  for  a  turbulent  or  transitional, 
boundary  layer  by  means  of  an  assumption  on  the  rate  of  mixing  with  th 
outer  flow.  In  their  work  the  boundary  layer  is  replaced  by  a  quasi -one 
dimensional  flow  with  the  aid  of  quantities  calculated  by  the  Prandtl 
boundary- layer  theory.  There  is  no  conflict  between  the  affirmative  con 
SEioS  of  the  present  paper  on  the  applicability  of  the  Prandtl  theory 
and  the  position  taken  in  the  Crocco-Lees  paper .J 
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9 

APPENDIX .A 

SYMBOLS 

The  following  symbols  are  used  in  this  report: 

a-l(x),  b1(x) 

Fourier  coefficients  appearing  in  the  solutions  for 
the  supersonic  stream 

e2Cx),  f2(x), 

G200,  H2(X) 

Fourier  coefficients  appearing  in  the  solutions  for 
the  subsonic  stream 

Designation  for  nine  integrals  considered  in  the 
analys is 

b 

/  thickness  of  subsonic  region 

f 

function  of  x  and  y  representing  incoming  waves 

g 

function  of  x  and  y  representing  outgoing  waves 

M 

Mach  number 

P 

static  pressure 

Ap 

change  in  static  pressure  ■ 

R 

reattachment  point  in  figure  5 

S  = 

RbXpgJ  also  separation  point  in  figure  5 

U 

undisturbed  velocity 

u,v 

disturbance  velocity  components  in  ,x-  and  y-directions 
respectively 

w  = 

x  +  p:Ly 

x,y 

rectangular  Cartesian  coordinates 

z  = 

x/r 

Pi  = 

-/%2  -  1 

p2  - 

V1  ■  m22 

r 

parameter  introduced  in  equation  (B3) 
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e 


A 

€  = 

!(«) 

P 

¥ 


deflection  angle  of  compression  wave 

displacement  of  interface  from  undisturbed  position 

Mo2 


density- 

disturbance  velocity  potential 


Subscripts : 
L 

U 


value  of  given  function  evaluated  at  the  lower  interface 

y  V-b 

value  of  given  function  evaluated  at  the  upper  inter¬ 
face,  y  =  0 


y  =  +o 

limiting  value  of  given  function  as  y  approaches 
zero  through  values  greater  than  zero 

o 

i 

ii 

limiting  value  of  given  function  as  y  approaches 
zero  through  values  less  than  zero 

l 

denotes  quantities  in  the  supersonic  region  y  >  0 

(excepting  Ip), 

2 

denotes  quantities  in  the  subsonic  region  -b  <  y 
(excepting  Ig) 

Superscripts : 

* 

indicates  limit  of  given  function  as  T  **■  0;  for 
example,  f(lf)*  =  li111 

r-*o 
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APPENDIX  B 


ANALYSIS 

In  the  present  report,  the  boundary  layer  has  been  replaced  by  a 
uniform  subsonic  stream  of  finite  width  b  which  is  bounded  on  one  side 
by  fluid  at  rest  and  on  the  other  side  by  a  uniform  supersonic  stream 
of  semi-infinite  extent.  For  purposes  of  analysis,  this  configuration 
will  be  considered  oriented  with  respect  to  a  rectangular  Cartesian 
coordinate  system  such  that  the  positive  x-axis  is  parallel  to,  and  in 
the  direction  of,  the  undisturbed  flow  and  coincides  with  the  undisturbed 
interface  between  the  subsonic  and  supersonic  streams,  and  such  that  the 
origin  denotes  the  point  of  incidence  of  the  given  shock  wave  with  this 
interface  (see  fig.  l) .  Thus,  for  -b  <y  <  0  the  flow  is  subsonic, 
while  for  y  >0  the  flow  is  supersonic.  Furthermore,  all  quantities 
in  the  supersonic  region  will  be  denoted  by  the  subscript  1  and  all 
quantities  in  the  subsonic  region  will  be  denoted  by  the  subscript  2. 


Supersonic  Region 

Governing  equation  and  solution  for  incoming  waves.  -  In  the  super¬ 
sonic  stream,  the  linearized  differential  equation  for  the  disturbance 
velocity  potential  <p-j_  is  given  by 


0 1 


d2^ 


dx‘ 


dy£ 


=  0 


(Bl) 


and  has  as  its  solution 

=  f(x  +  pxy)  +  g(x  -  pxy) 


(B2) 


where  f  and  g  are  arbitrary  functions.  The  undisturbed  flow  has 
been  taken  in  the  positive  x-direction  so  that  waves  represented  by 


f(x  +  p^y)  will  be  incoming  waves  while  those  represented  by 
g(x  -  p^y)  will  be  outgoing  waves  emanating  from  the  upper  i; 


For  purposes  of  the  analysis,  the  form  of  the  function  f 
specialized  and  taken  to  be 


interface . 
will  be 


f(x  +  Ppy)  = 


-r(x  +  pxy) 


4 


x  +  Ppy  <  o 


x  +  P2y;  >  0 


(B3) 


where 

•where 


U-j_  is  the  undisturbed  velocity  in  the  supersonic  stream  and 
X  and  e  are  parameters.  If  y  tends  to'  zero,  the  incoming 


waves  will  degenerate  into  a  single  compression  wave  with  deflection 
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angle  e,  and  the  point  of  incidence -of  this  single  compressive  wave 
with  the  upper  interface  will  he  the  origin  of  the  given  rectangular 
Cartesian  coordinate  system.  - 

Now,  by  definition,  the  components  u  and  v  of  the  disturbance 
velocity  are  given  in  terms  of  the  disturbance  velocity  potential 
by  the  relations 

(B4) 

oy 


_  ^3. 

3x* 


Therefore,  the  components  of  the  disturbance  velocity  of  the  incoming 
waves  are  given  by 

r  ' 

0  ,  •  -x  +  0xy  <  0 


ul, incoming  < 


(B5) 


-Uxe  -y(x  +  p^y),  x  +  pxy  >  0 

^~0Te 


and 


r 


'l, incoming 


0 


-U-jee 


-T(x  +  pxy) 


X  +  P-^  <  0 


x  +  p^y  >  o 


(B6) 


'  '  Incoming  and  outgoing  waves  in  Fourier-integral  form.  -  It  will 
prove  useful  to  express  the  disturbance  velocities  of  the  incoming  waves 
in  the  equivalent  Fourier- integral  form.  To  do  this,  it  is  noted  that 
the  components  of  the  disturbance  velocity  of  the  incoming  waves  are  ^ 
proportional  to  the  function 

(o  ,  w  <  0 


h(w}=\ 


(B7) 


vf 


-yw 


w  >  0 


where  v  =  x  +  pxy,  and  that  h(w)  can  be  represented  in  an  entirely 
equivalent  form  by 


cos  Xw  + 


— — tt  sin  Xw)  dX  (B8) 


^  .2 
y  +  X 
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Furthermore,  the  potential  due  to  the  outgoing  waves  can  also  he  expressed 
in  Fourier -integral  fom,  namely, 

g(x  -  pxy)  =  Ux  j  [^(X)  sin  X(x  -  p^y)  +  B^X)  cos  x(x  -  Ppyjj  dX 

(B9) 

where  Ap(X)  and  B1(x)  are  undetermined  Fourier  coefficients.  Thus, 
after  differentiation  of  the  expression  for  g(x  -  p^)  to  obtain  the 

disturbance  velocities  due  to  the  outgoing  waves,  the  separate  distur¬ 
bance  velocities  can  be  superposed  and  the  Fourier-integral  expressions 
for  the  components  of  the  total  disturbance  velocities  (of  both  incoming 
and  outgoing  waves)  obtained.  Evaluated  at  the  upper  interface,  these 
disturbance  velocity  components  will  then  be  given  by  ^ 


l,y=+0  irp-L 


— _  cos  +  -j 

r2  +  X2  r2 


— — 5-  sin  Xx  )  dX  + 
+  x  / 


‘/t 


^1  /  |XAj_(X)  cos  Xx  -  XB1(x)  sin  Xx  dX 


(BIO) 


%y=+0  =  HT  J  (“T cos  Xx  +  ^77  sin  Xx) 
PA  (X)  cos  Xx  -  XB1(X)  sin\xj  dX 


cos  Xx  + 


(Bll) 


Subsonic  Region 

Governing  equation  and  general  solution  in  Fourier-integral  form.  - 
For  the  subsonic  stream,  -b  <  y  <  0,  the  linearized  differential  equation 
for  the  disturbance  velocity  potential  is  given  by 


2  d2<P?  S2cp, 


(B12  ) 
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and  the  solution  <p2  can  Be  expressed,  among  many  possible  forms,  as 
<p2(x,y)  =  U2^/^  cosh |x.j32(y+b^j  *  |e2(X)  sin  Xx  +  F2(X)  cos  Xxj  dX  H 


sinh  |^|32(y+b)|  '  j^2 (^)  sin  Xx  +  H2(X) 


cos  Xx  dX 


(B13) 


where  E2(X),  F2(X),  Gg(X),  and  H2(X)  are  also  undetermined  Fourier 
coefficients.  Then,  by  retracing  the  steps  just  taken  for  the  super¬ 
sonic  case,  there  is  obtained 


u2(x,y)  =  U2  I  coshjxp2(y+b")j  •  |e2(X)  cos  Xx  -  F2(X)  sin  Xxj 


X  dX  + 


sinh jxp2(y+b)j  •  ^2 C^-)  cos  ^  “  H2(X)  sin  Xxjx  dX 

(B14) 


v2(x,y)  =:  U2p2  I  sinh|xp2(y+b)j  •  |e2(X)  sin  Xx  +  F2(X) 


cos  Xx  X  dX  + 


U2p2  /  coshjxp2(y+b)J  *  jc2(X)  sin  Xx  +  H2(X)  cos  Xxjx  dX 

Jo  (B15) 


Boundary  Conditions 

In  these  solutions  for  the  subsonic  and  supersonic  streams,  six 
unknown  Fourier  coefficients  Ai(X),  BX(X),  E2(X),  F.2(x),  G2(X),  and 
22(X)  exist.  However,  these  coefficients  are  determinable  by  boundary 
conditions  at  the  upper  and  lower  interfaces.  The  physical  conditions 
which  must  be  satisfied  across  these  interfaces  are:  (l)  equal  static 
pressure,  and  (2)  equal  inclination  of  the  flow.  At  the  upper  interface, 
condition  (l)  gives,  in  the  small  perturbation  approximation, 
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(I)  =  -P2U2^2,y=-0  for  a11 


x 


.  (B16) 


At  the  lower  interface,  the  perturbation  pressure  is  taken  to  he  zero; 
condition  (l)  therefore  reduces  to 


(II)  u2>y=_b  =  0  for  all 


(B17) 


Condition  (2)  gives  no  information  at  the  lower  interface  since  the  flow 
speed  is  zero  below  the  interface;  at  the  upper  interface,-  this  condition 
yields  (small  perturbation  approximation) 


(III)  v2>y=s_o  Ux  =  v]_ } y=+0  U2  for  arL  X 


(B18) 


Establishment  of  the  four  equations  for  the  four  Fourier  coefficients 
With  use  of  equation  (B14),  boundary  condition  II  implies 


cos  Xx  -  F2(\)  sin  Xxjx  dX  =  0 


(B19) 


so  that  E2(X)  and  F2(X)  vanish  identically  and  <p2  reduces  to 

sinh  |xp2(y+b)]  [02(X)  sin  Xx  +  H2(X)  cos  XxJ  dX 


(B20) 


u2  =  U2 


to 


sinh  ^p2(y+l3)]  [>-G2 Cx )  cos  Xx  -  XH2(x)  sin  Xxj  dX 


(B21) 


and 


v2  =  U2P2  /  coshjxp2(y+b)j  |xg2(X)  sin  Xx  +  XH2(X)  cos  Xxj  dX 

(B22) 

Substituting  equations  (BIO)  and  (B2l),  the  latter  evaluated  at 
y  =  -  0,  into  equation  (B16 ) ,  that  is,  boundary  condition  I,  gives 
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Piui 


(yb?  “  yh?  “  xj  awJ  [x*ltk)  “■* '  “lW  ^ ll) 

=  P2U ^  f  sinh  (X.bp2)|xG2(X)  cos  **  ‘  XH2(X)  sin  Xxj  dX  ,  /  . 


(B23) 


whence,  upon  equating  corresponding  Fourier  coefficients  and  making  use 
of  the  fact  that 


ppUo2  M?2 

_22  =  2  for  =  p2 

PlUi2  MX2 


there  is  obtained 


Ai(\)  -  G2(X)  “2  sirJh  (^2)  =  ^  r2  l  x2  X 


and 


M2 


2 


Bi(X)  -  H2(X)  ^  sinh  (Xbp2)  =  %.-J7^2 


-e_  1 

r 


Similarly,  boundary  condition  III  (equation  (B18)')  implies 


(B24) 


(B25) 


(B26) 


U1U2P2  1  COSh 
' 0 


(\bp2)  jxGg(X)  slrl  **  +  MzW  003  dX 


=  *Mf)l  hf7C0B^  +  rz 


b—-g  sin  Xx^  dX  -  |XAX(X)  cos  Xx  -  XBX(X)  sin  Xx]  dxj 


so  that 


AX(X)  +  H2(X),  ^  cosh  (Xbp2)  =  -2^2 


(B27  ) 
(B28) 
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Bl(X)  -  G2(X)  pj  cosh  (Xbp2)  -  ^  +  ^2 


(B29  ) 


Hence ,  the  prohlem  is  now  reduced  to  one  of  evaluating  the  four 
undetermined  Fourier  coefficients  Ai(\),  B1(\) ,  G2(X),  and  Hg(X)  by 

means  of  the  four  equations  (B25),  (B26),  (B28),  and  (B29) . 


Evaluation  of  Fourier  Coefficients 

Solving  equations  (B25),  (B26),  (ESS),  and  (B29)  simultaneously 
yields  the  following  expressions  for  G2(X)  and  H2(k): 


"/m  2\2  2  "I 

(-^5-)  sinh2(Xbpo)  +  cosh2  (Xbp2) 

IWJ  h2  J 


+  |p  cosh  (Xbp2) 


(B30) 


H2(X)  cosh2  (Xbp2)  + 

Jl 


sinh2  (Xop2)| 


*Pl  ‘  r2 


1  p2 

+  X2 


Y  cosh  (Xbp2) - o  sinh  (Xbp2) 

a.  “  Mi*- 


(B31) 


In  addition,  if  the  Mach  number  parameter  9  is  defined  such  that8, 

M22 


9 

cos  ^  = 


422\2  P22 

2  )  +  2 
i  \k  &  J  n 


^Although  simplicity  is  not  necessarily  thereby  achieved,  0  is  intro¬ 
duced  so  as  to  maintain  a  parallel  presentation  with  that  of  reference  2 
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then 


©’ 


p22  1  /M22\  cosh(2Xbp2)  -  cos  9 

inh2  (Xb02)  +  ^2  COsil  CXbP2)  =  2  cos2'£ 


Pi 


(B33) 


so  that  G2(X)  and  H2(X)  are  also  given  hy 

2 

h  /Mo^ 

g2(x> 


i  /V  \2 

_2  W2/ 


cosh  (2Xb02)  -  cos  0 


2  ® 
COS  2 


I 


-2  e  1 

*Pi  r2  + 


r  /m22\  . 

S1 


P2 


inh  (Xhp2)  +  —  cosh  (Xb02) 


■3 


(B34) 


and 


1  cos^  (2Xbp2)  -  cos  9 

B-zW  1 


vMl2/ 


cos 


2  £ 


-2 

«0 


e  i  l"P2  r 

i  'r2  t  x2  Ifi  x 


m2£ 


cosh  (Xb0~)  - — ^  sinh  (Xb02) 
c  M-/ 


>5 


(B35) 


Similarly,  hy  means  of  the  same  equations,  the  following' expressions  for 
A-JX)  and  Bj^X)  are  also  obtained: 


k±{\)  cosh2  (Xhp2)  +  (%^j  Sinh2  (Xbfig 

,  2  ^  "1 
=  -I-  5  1  g  K  Y  cosh2  (Xbf»2)  -  (^)  l  sinh2  (Xhp2)  -  2  ^2  t  sinh  (Wj'1  cosh 

j%  r2  +  xZ  l_Pl2  X  •  \Ml  /  %  pi  J 


(B36) 
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and 


B1(X)  sinil2  (Xbl32),+  ~~2  cosh2  (^2^ 


—  #  1  r  6p  o  Bp^  o 

^p7  "2  —  -jj-g  £  X  j^2  pj  sinil  CX‘bp2)  cosh  (Xbp2)  -  sinh  +  7T  cosh2  (XbP2) 


(B37) 


Solution  for  Pressure  Along  Upper  Interface 


With  the  Fourier  coefficients  Ai(X),  Bl(X),  G2(X),  and  H2(X) 

evaluated,  the  solution  of  the  originally  stated  problem  is  now  complete 
However, ' there  exist  other  quantities  of  interest  such  as  the  pressure 
distribution  along  the  upper  interface  and  also  the  displacement  of  each 
interface  from  its  undisturbed  position.  These  quantities  will  now  be 
considered. 


If  the  change  in  pressure  from  the  undisturbed  value  is  denoted  by 
Ap,  then  the  linearized  value  is 


Ap  =  -  pUu 


Hence,  denoting  quantities  at  the  upper  interface  by  the  subscript  U 
and  those  at  the  lower  interface  by  the  subscript  L,  equation  (B2l)  and 
-  boundary  condition  I  give 


(AP)y  =  -P2u2u2,y=-0  =  ~P2U2  2  sinh  (Xbjig)  ^XGg(x)  cosXx  -XHg(X)  sin  Xxj  dX  (B38) 


where  G2(\)  and  ^(x)  are  given  by  equations  (B34)  and  (B35).  Or, 
with  use  of  equation  (B24)  and  also  since  equation  (B32)  further  implies 
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0  _ 
sm  g  = 


rl  -  cos*  £ 


>2\2  P22 


(B40) 


equation  (B38)  can  then  he  written  as 


to)u  4*  f  2  8 

1  ti2  *Pi  .  .  2 

2  plUl  L  .  i 


(777  cos  Xx  +  ?-,?  '':r'  ^  dX  ’ 


I  sln2  9 


y  cos  Xx  +  X  sin  Ax 
(r2  +  x2)(cosh  2Xbpg  -  cos  0) 


I sln  9 


X  sinh  (2XbP2)  cos  Xx  -  y-  slnh  C2XbP2)  sln  ** 


(y2  +  \2)(cosh  2Xbp2  -  cos  0) 


(B41) 


Limit  as  r-»0  (incoming  waves  -» step  function).  -  If  T-*-0,  the 
incoming  waves  will  degenerate  into  a  single  compress ionwave  with 
deflection  angle  e.  Consequently,  in  order  to  obtain  the  value  of  the 
pressure  change  at  the  upper  interface  for  a  single  compression  wave  of 

(Ap)u 

deflection  angle  t,  let  r-0  in  the  expression  for  T 

2  plUl 

From  equation  (B8)a 


^ere,  the  subscripts  under  I  designate  the  various  integrals  rather 
than  any  position  in  a  subsonic  or  supersonic  stream. 
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!l  = 


'0 


( 


,  x  <  0 

- — cos  Xx  +  — — — -  sin  Xx]  dX  =  ^  %  ,  x  =  0 


r2+*.2 


2  -.2 

r  +  x 


iin  Xx^ 


rte^*,  x  >  0 


(M2) 


whence 


cos  Xx  +  —7^— — 5-  sin  Xx)  dX  =  <  -|  , 


r2  +  X2 


■) 


0  , 


x  <0 
x  =  0 
x  >  0 


(M3) 


For  the  other  integrals  of  equation  (B4l),  a  theorem  justifying  the 
interchange  of  the  processes  of  integration  and  of  passing  to  a  limit  is 
necessary.  One  such  appropriate  theorem  which  will  suffice  for  the 
analysis  in  this  report  will  now  be  stated  (reference  10,  p.  474). 


DEFINITION .  -  Let  R  .=  (a  «>  a# )  denote  the  rectangle  determined  by 
a  <  x  <  «,  a  <  t  <  6,  3  finite  or  infinite.  The  function  f(x,t)  is 
said  to  be  REGULAR  IN  R  when:  (l)  f(x,t)  has  no  point  of  infinite 
discontinuity  in  R;  (2)  f(x,t)  is  integrable  in  a  <  x  <  »  for  each 
t  in  a  <  t  <  p . 

THEOREM  I.  -  (l)  Let  f(x,t)  be  regular  in  R 
or  infinite,  except  possibly  on  the  lines  x  =  . 

f(x,t)dx 

be  uniformly  convergent  in  a  <  t  <  p.  (3)  As  t-+tQ,  tQ  finite  or 

infinite,  let  f(x,t)  approach  a  limiting  function  F(x)  uniformly  in 
any  region  a  <  x  <  b,  except  possibly  on  the  lines  x  =  a^  .  .  . ,  x  =  ar 

Then,  p» 

lim  /  -  f(x,t)  dx  exists. 
t-t0  / 


=  (a*»ap),  p  finite 
.  .,  x  =  ar.  (2)  Let 
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(4)  Let  F(x)  be  integrable  in  any  region  a  <  x  <  b, 


Then, 


lim  /  f (x,t)dx  =  /  F(x)dx 

t  -+-tn  / 


Thus,  for 


_ y  cos  Xx 

(y2  +  X2)(cosh  2bXpg  -  cos  ©) 


(B44) 


first  letting  z  =  X/V  so  that 


cosiyxzj 


(z2  +  l)(cosh  2byzP2  -  cos  0 ) 


(B45) 


and  then  applying  theorem  I,  there  is  obtained 


lim„  x2  =  1  -  cos~e' 
y-*  o 


z2  +  !  2(1  -  COS  e4 


4  sin' 


•  2  y 
i  rr- 


(B46) 


*3  f 


X  sin  Xx 


Q  (r2  +  X2) (cosh  2bp2X  -  cos  9) 


(B47) 


X  sinh(2bp2^)  cos  Xx 


(y2  +  x2)(cosh  2b0-2X.-  cos  0) 


(B48) 
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let 


s  —  25X62 


K  = 


rtx 

2b  p,. 


so  that 


G  •  Sg 
S  sin  — 


fr2  + 


4b2p22 


) 


dS 


(cosh  S  -  cos  0) 


and 


S  sinh  S  cos 


I4  = 


Jt 


dS 


[r2  +  — ! — o  ]  (cosh  S  -  cos  0) 

\  4b2(32  / 


Then, 


sin 


SZ 


S(cosh  S  -  cos  0) 


dS 


and 


lim  I4  = 
r-*-  0 


S  E 

sinh  S  cos  — 

_ it 

S(cosh  S  -  cos  0) 


dS 


'0 


Lastly, 


X  sinh(2Xbp2)sin  Xx 
(y2  +  \2)(cosh  2Xb(32  -  cos  0) 


d\  = 


Solution  for  single  incoming  compressive  wave.  -  Combining 
five  limits • yields 


23 

(B49) 

(B50) 

(B5l) 

(B52) 

(B53) 

(B54  ) 

0 

(B55) 

these 
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(Ap) 


U 


ilm  r  2 

r-°  z  Piui 


4 e 

Pi 


cos2  | .  l(C )  - 


|  sin20 


\4  sin^  of 


Bin 


SK 


*?,  n  I 

2jt  Sin  6  /  S^cosh  s'  -  cos  0 ) 
>0 


where  l(€)  Is  defined  by 


r 


o, 

i(C) 


V 


1> 


£  <  o 
£  =  0 
?  >0 


(B57) 


Now,  it  can  he  shown  by  contour  integration  that,  for  £  >  0, 


.in  f  __  l  _  , 

2it  /  S(cosh  S  -  cos  0)  ^  i  j  fl  +  2«  sin  0  ^ 

*  sin  o 


£ 

n=0 


.  -(2n+l)E 


0  6 
2n  +  2  -  ~  2n  +  j- 


?>0 


(B58) 


and 


sinh  S  cos  -y-  1  -(2n+l)£ 

Jt  dC=  —  /  e 

2it  /  S(cosh  S  -  cos  0)  2it  ^ 

/0 


-fir> 

+ 


0 


0 


2n  +  2  -  „  2n  +  ^ 


,£  >0 


J(B59) 


Furthermore,  since  the  integrands  are  odd  and  even  functions  of  £, , 
respectively,  it  can  immediately  be  deduced  that,  for  £  <  0, 


sinf  -i 

dS  =  - 


2n  ,  '/  S(cosh  S  -  cos  ej  8  sinn.  x. 


,2  £ 

2 


.  1.-V 

2«  sin  9  /  . 


(2n+l)? 


n=0 


■  ^M_s 


2n  +  2  -  -  2n  +  £ 
«  « 


,  e<o 


(E60) 
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and 


sirih  S  cos  -5- 

_ it 

S(cosh  S  -  cos  0) 


dS 


2*  t—t 
n=0 


(2n+l}5 


2n  +  2  - 


2n  +  - 
Jt 


j  5<0  (B6l) 


Therefore,  with  the  substitution  of  these  results  into  equation  (B56), 
it  follows  that: 


r 


(Ap)tj 

lim  1  2 
l"*-0  2  PlUl 


4c  .  . 

-r—  sin  0 

*P1  5^6  2n  +  2 


/  /  q  >  K  <  o 

1  1.  o  n  w 


4e 

«Pl 


sin  0 


•O 

E 

'  n=0 


2n+2 


2n  +  — 
Jt 


,  C  >  0 


(B62) 


(B63) 


Displacement  of  Upper  Interface 

A  quantity  of  primary  interest  is  the  shape  of  the  distorted  inter¬ 
face  between  the  supersonic  and  the  subsonic  flows .  The  slope  of  this 
interface  is  related  very  simply  to  the  already  determined  pressure  at 
the  interface,  according  to  the  following  considerations.  If  the  dis¬ 
placement  from  the  undisturbed  position  is  denoted  by  r)|j,  this  slope  is 
given  by  - 


drju  ^  V2 
dx  Uo 


all  £,  y  =  0 


(B73) 


approximately.  Now,  applying  the  boundary  condition  of  equal  slopes  on 
both  sides  of  the  interface,  equation  (B18)  gives 


Vn 


Uo  Un 


all  z,  y 


(B74) 


Furthermore,  v^  can  be  related  to  uj_.-  In  the  region  £  <  0,  only  out¬ 
going  waves  exist,  corresponding  to  the  g-function  in  equation  (B2);  thus 
by  differentiation  of  g,  , 
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V1  =  "  HU1  C  <  0,  y  =  0  (B75a) 

In  the  region  £  >0,  there  exists,  in  addition,  the  incoming  weak  shock 
providing  a  flow  deflection  e  for  r-*0,  corresponding  to  the  f -function 
in  equation  (B2);  this  provides  an  additional  term  such  that 

lim  vx  =  -  (3^  -  2eU1  £  >  0,  y  =  0  (B75h) 

T*  0 

When  these  equations  are  combined,  there  results 


lim  — r— 
dx 

Y-+  0 


ui 

‘  P1  U-i 


K  <  0,  y  =  0 


(B76) 


-  Pi  ir  " 2e,  K  >  °>  y  =  0 


Then,  by  virtue  of  the  relation  Apy  =  -  PpUiui^  this  may  be  written 


y*  0 


1  Q  I  APU 

2  Pl  I.  .2 

\2  plUl  , 


K  <  o,  y  =  o 


2  H 


(B77 ) 


2«,  £  >  0,  y  =  0 


This  is  the  desired  relation  between  the  local  slope  and  the  pressure  of 
the  interface., 

Upon  substitution  of  the  calculated  pressure,  equations  (B62)  and 
(B63),  there  results 


-  (2n+2  - 

V*  ej _ *]_ 

L/  ,  e 


2e  \  e ' 

—  sin.0  7T~~e 

n=0  ,  2n  +  « 


g  <  0 


(B78) 


,  .  <% 

lam  — — 

.  dx 
y*  0 


•  9  V  e'C2n  +  0e 

j x  Lj  ~  .  e 

n=0  2n  +  7 


5  >  0 


(B79) 
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Expressions  for  displacement  of  upper  interface.  ^Denoting 


d%  ,  *i„* 

by  -ST 
r-*-o 


,  equation  (B50)  gives 

dV*_  dV*d£_  d%*  n 

dx  d£  dx  d£  2bp2 


(B80) 


Thus,  with  the  assumption  that 


%*{-•;  -  o 


(B8l) 


dTJu* 

the  expressions  for  ■  can  be  integrated  and  will  yield 


^n+2  -  ^ 

%*fei  9)  -%■  sin  e  Y,  J- - -T2->  «  <  0 

n'°H-|)  (B82 ) 


while  for  £  >  0,  since 


<%*(?<  0)  0)  /  drlu*(?>0;  0) 


d^ -  d?  (B83) 


there  is  obtained 


V*te>  0)  =  sin  0  X) 


n=0  fa  ~  |) 


i _ -  -  2e  sln  0  V  £ 

A  2  In  . 


-  -(2n+D« 


(2n  - !)‘ 


(B84  ) 


Hence 


•  -fzn  +  £)l 

Vfe*  "  T  Sl“  9  X/  "7 - - - 2  + -  1  g  -  5 - 2 

-=0  (2„2-§J  (2n  +  |)  (&.  +  §) 


-  2e«, 


>  „  <Be5> 
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Therefore,  in  summary, 


ZbB„  0) 


2e 

m 

V  , 

r*4  m  a  / 

* 

sin  C7  /  j 

n=0 

(2n+2  - 1)2 

«• 

2e 

i 

JC 

sin  u  /  / 

(2n+2  -  f 

1*4 


£  <  0  (B82) 


2e£,  £>°  (B85) 


By  inspection  of  equations  (B82)  and  (B85),  it  is  noted  that 
%*(£}  9 )  is  finite  and  continuous  at  '  £  =  0,  while,  on  the  other  hand, 
equations  (B78)  and  (B79)  show  that  di^/dx  is  infinite  at  %  =  0. 


Displacement  of  Lower  Interface 

Finally,  the  expression  for  the  displacement  of  the  lower  interface 
will  he  deriyed.  If  the  displacement  from  the  undisturbed  position  is 
denoted  by  t]l,  this  displacement  will  have  a  slope  given  in  linearized 

approximation  by  . ■ 


d^L  v2,y=-b 
dx  Uo 


(B86) 


or,  with  use  of  equation  (B22), 

dlT  r*  r  . 


TE-/  W«.{0*tXj  sin  Xx  +  H2(X)  cos  Xxj  dX 


(B87  ) 


I 

where  G2(X)  and  H2(X)  are  given  by  equations  (B34)  and  (B35).  That  is, 
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Limit  as  r4  0  (incoming  wave  -» step  function) .  -  Consideration  will 
now  be  given  the  various  parts  of  the  expression  for  dTjy/dx  as  one  again 

passes  to  the  limit  y  =  0.  For 


I6  5 


y  sinh  (Xbpo)  sin  Xx 
j  - : - dX 

Jq  (y2  +  \2)(eosh  2bX^2  -  cos  0 ) 


(B89) 


first  letting  X  =  yz,  there  is  obtained 

sinh  (ybzpg)  •  sin  (yxz) 


lim  Ig  =  lim  . 

-p*  0  y-*.o  /  (l  +  z^)(cosh  2yzbp2  -  cos  9 ) 

f0 


dz  =  0  (B90) 


Similarly,  for 


Ir,  s 


y  cosh  (XbPg)  cos  Xx 
(y2  +  \2)(cosh  2Xbpg  -  cos  9) 


dX  (B91) 


lim  I- 
y-*-0 


cosh  (yzbpo)  cos  (yzx) 
lim  /  dz 

y-0  /  (1  +  z2)(cosh  2yzbP2_-  cos  9) 

to 


(B92) 


1  /  dz 

to 


_  1  rt 

1  -  cos  9  I  ,  ,  _2  ’  1  -  cos  0  2 

1  +  z 


/ 


30 


so  that 


NACA 


lim  X7  =  ,  . . ^0 

0  -  .4  sm  g- 


For 


Is  '  '  2  ,  \  2 


X  cosh  (Xb(32)  •  sin 


ax 


and 


X9S 


(y2  +  X2)(<?osh  2Xbj32  -  cos  0) 


X  sinh  (XToPg )  .  cos  Xx 
(y2  +  X2)(cosh  2Xbpg  -  cos  0) 


dX 


introducing  S  and  £  as  before  yields 


*8  = 


n  v.  S  •  §£ 

S  cosh  o  sm  -f 
^  a 


arid 


/o  (^2  +  4b2  t)  (cosh  S  -  COS  0) 
^  2 


-  .  ,  S  S£ 

s  sinh  ^  cos  -f 

b  7t 


dS 


%  " 


2  \ 

/0  (r2  +  ” i — 2)  (cosh  s  -  cos  0) 

\  4b  Bo  ’ 


dS 


4b2p22 


whence 


lim  In  = 
y-+  0 


cosh  75-  sin  — ^ 

- 2 - jt -  ds 


S  /  S(cosh  S  -  cos  0) 


and 
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(B93) 


(B94) 


(B95) 


(B96) 


(B97  ) 


(B98) 


(B99) 
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Furthermore,  following  reference  2,  it  can  be  shown  that,  for  K  >  0, 


r 


.  .20  «  .  V  t t 

4  sin  j  2  sm  ^  n=o 


lUn  I8  =  < 
r->o 


4  sin2  |  2  sin  |  feo 


^  (-I)"  e-(2n+l)5 


„  „  e  „  e 

2n  +  2  -  —  2n  +  - 

rt  jt 


:-ra«  ,e# 


2n  +  2  -  —  2n  +  — 
jt  jt 


K<  0  (B100) 


£>0  (B101) 


and. 


1  £  (-I)"  *(2n+l)5 


2  cos  —  n=0 


.-69* 


lim  I  =  < 
r-»o  s 


i-jE  w"  «'t2”1|E 


2n  +  2  -  —  2n  +  — 
rt  rt 


2  cos  -  n=0 


V. 


2n  +  2  -  -  2n  +  - 

jt  Jt 


,  K<o 


>  K>o 


(B102) 


(B103) 


Simplification  of  expressions  for  slope.  -  Combining  these  limits 
and  using  equations  (B32)  and  (B40)  yield 


lim 

p -+  0 


dTlL 

dx 


2e  sin2  | 


2  ^sin2  |-  Jt  sin  n=0 


V  -(2n  + 1) 

L,  C-h  e 


2n  +  — 
jt 


b  K>0 


(B104) 


That  is. 


lim 

If  -*■0 


d\ 

dx 


-  2e  1 1 - sin 

Ti  O 


*  -4n  +  A 


b  K  >  o 


2n  + 


it 


(B105) 
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Furthermore ,  by  symmetry  considerations  of  the  various  integrands,  it 
follows  that,  for  £  <  0, 


dTir 

lim  — =  -  2e  sin 
_  dx 

X  -►O 


2  1  jt  sin  ^  n=0 


+. 


2n+2  - 


it 


2n  +  - 

It 


(-I)”  e(2n+l)5 


it  sin  g  n=0 


e(^)«  e-M? 


2n+2  -  £  2n  +  — 


That  is. 


drir 

lim  ~  =  -  2e 

T-+0 


I sin  l  S  ^'i^n 

*  2  feo 


(2n+2  - 


2n  +  2  -  - 
It 


,  g  <  0  (B107) 


Therefore,  in  summary, 

r  f 

'  /2  -  6 
-2e  I—  sin  o 
Art  *  £ 


E  (-r 

n=0 


k 


n+2  -  — )? 


e  1  ’ 

2n  +  2  -  - 

it 


dlL  y 

~  ^  /  ■  A  -L  ,  2> 

L  2  .  .9  V  ,  ,  ,n  e  V  V 

I1  *  ?  S“  S  U  C'1).  2n  +  £  '  ’ 


y-*-0 


-26  I 


£<  0  (B107  ) 


£  >  0  (B105) 


n=0 


Final  expressions  for  displacement  of  lower  interface.  -  If 

dr)  dTW* 

lim  -g=  is  denoted  by  -r— ,  it  is  noted  that 

r-o 

d\*  d\  3T  aV  £ 

dx  d£  dx  d£  ’  2bp2 


(B108) 


so  that  integration  of  equations  (B105)  and  (B107),  along  with  the  boundary 
condition  > 


%*(-»?,  0)  =  0 


(B109) 
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yields,  for  %  <  0, 


2bp 2 


r/zn+2  - 

c‘1)n  r~7\ 

H-!) 


2~  )  K  <0 


(B110) 


For  £  >  0, 


-0 


T)L*(Si0)  = 


dr]L*fe<O;0) 


d£  + 


di]T*(5>O;0)  . 


d? 


so  that 


ni*(C;9)  =  -ae  +  £  Bln  |  J  (-l)n 


n=0 


-1 


1  e 

rs+  ~ — ^ 


rh  *  D' 


M-!)  M)  (2“  ♦  I) 


,  K>o 


(Bill) 


Therefore,  in  summary 


P2 


a  H  - 

-  t  sln  1 1/  C'1)11  - - —  «  «  <  0 


n=0 


h*  ■  I) 


Zbpi  =  S 


-■UK 


4E  .  0  )  /  ,  *n 

tTsul2^  (_1) 


n=0 


.  -1 


*  §« 


(B110) 


,  e>o  (Bin) 


Both  t}l*(£;0)  and  dr)L*/dx  are  continuous  at  5  =0.  The  con¬ 
tinuity  of  t)l*{  5; 0 )  is  readily  seen  while  the  continuity  of  dr^/dx 
may  he  shown  by  letting  x  =  0/ 2*  in  the  formula  for  jt/sin  Jtx  given  on 
page  208  of  reference  10.  * 
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TABLE  I  -  PERTURBATION  OF  UPPER  INTERFACE 


*  V 

zb^~n  against  •« 


7TX  . 

2bp2 


for  6  =  rt/4,  rt/2,  and  3rt/4 


4 

&=rt/4 

9=itf2 

a=3rt/4 

-5.0 

0.00001 

0.00008 

0.00028 

-4.0 

.00007 

.00035 

.00097 

-2.5 

.00093. 

. 00335 

.00632 

-1.75 

.00345 

.0103 

.01620 

-1.00 

.0132  . 

.0324 

.0416 

-  .75 

. .0208 

.0480 

.0584 

-  .50 

.0335 

.0722 

.0820 

-  .40 

.0410' 

.0856  ■ 

.0944 

-  .25 

.0562 

.1120 

.  1180 

-  .10 

.0804 . 

.  1520 

.1517 

-  .05 

.0925 

.1715 

.  1675 

0 

.1108 

.1993 

.1891 

.01 

.1138 

.2009 

.1856 

.05 

.1215 

.2018 

.1679 

.10 

.1253 

.1958 

,1408 

.25 

.1240 

.1564  •. 

.0450 

.40 

.1169 

.1003 

-0.0620 

.do 

.1017 

.0566 

-  ,1377 

.75 

.0516 

-0.0686 

-  .3363 

1.00 

-0.0137 

-  .2117 

-  .5450 

1.75 

-  .2803 

-  .7168 

-1.2120 

2.50 

-  .6323 

-1.3002 

-1.9155 

4.0 

-1.5293 

-2.6075 

-3,3750 

5.0 

-2.2368 

-3.540 

-4 . 3650 
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TABLE  II  -  PERTURBATION  OF  LOWER  INTERFACE 


2bp2  2e 


-JEL-  for  9  =  jt/4,  rt/2 ,  and  3it/4 

2'd  P  2 


K 

9=rt/4 

-5.0 

-0.000012 

-4.0 

-  .00007 

-2.5 

-  .00100 

-1.75 

-  .00368 

-1.00 

-  .0134 

-  .75 

-  .0205 

-  .50 

-  .0309 

-  .40 

-  .0362 

-  .  25 

-  .0459 

-  .10 

-  *0577 

-  .05 

-  .0622 

0 

-  .0670 

.05 

-  .0721 

.10 

-  .0775 

.25 

-  .0959 

.40 

-  .1173 

.50 

-  .1338 

.75 

-  .1811 

1.00 

-  .2389 

1.75 

-  . 4735 

2.5 

-  . 7942 

4.0 

-1.6420 

5.0 

-2.3253 

S=rt/ 2 


-0.00011 

-  .00050 

-  .00472 

-  .0144 

-  .0436 

-  .0625 

-  ’.0890 

-  .1021 

-  .1250 

-  .1522 

-  .1623 

-  .1728 

-  .1841 

-  .1962 

-  .2358 

-  .2803 

-  .3130 
.4065 

-  .5168 

-  .9293 
-1.4453 
-2.6725 
-3.5775 


&*3it/4 


-0.00073 

-  .00254 

-  .01653 

-  .0420 
-' .1057 

-  11428 

-  .1918. 

-  .2150 

-  .2550 

-  .3005 

-  .3173 

-  .3345 
r-  .3528 

-  .2718 

-  .4330 

-  .5013 

-  .5503 

-  .6860 
-  .8385 
-1.3803 
-2.0100 
-3 . 4025 
-4.4675 


2665 


2665 


4.0 


(c)  Mach  number  parameter  9,  3jt/4. 

Figure  3.  -  Concluded.  Pressure  increment  Ap^  plotted  against  distance  x  from  point 

of  incidence  of  wave  upon  upper  interface. 
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Figure  4.  -  Displacement  of  upper  and  lower  interfaces  plotted  against  distance  from  point  of  incidence  of  wave  upon  upper  interface. 
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